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ABSTRACT 

Several X-ray spectral models for tori in active galactic nuclei (AGNs) are available to 
constrain the properties of tori; however, the accuracy of these models has not been 
verified. We recently constr uct a code for the torus using Geant4, which can easily 
handl e different geo metries dLiu k hill 20 Mil . Thus, we adopt the sa me assu mptions as 
iMurphv k Yaaoobl (I2009L hereafter MY09) and lBrightman k Nandral (1201 ll . hereafter 
BN11) and try to reproduce their spectra. As a result, we can reproduce well the 
reflection spectra and the strength of the Fe Ka line of MY09, for both TVh = 10 24 
and 10 25 cm -2 . However, we cannot produce the strong reflection component of BN11 
in the low-energy band. The origin of this component is the reflection from the visible 
inner wall of the torus, and it should be very weak in the edge-on directions under the 
geometry of BN11. Therefore, the behaviour of the reflection spectra in BN11 is not 
consistent with their geometry. The strength of the Fe K a line of BN11 is also different 
from our results and the analytical result in the optically thin case. The limitation of 
the spectral model will bias the parameters from X-ray spectral fitting. 

Key words: radiative transfer - galaxies: active - X-rays: galaxies. 


1 INTRODUCTION 

Under the unification scheme of active galactic nuclei 
(AGNs), a toroidal structure referred as “the torus” pro¬ 
vides anisotropic obscuration and can explain the di versit y 
of optical and X-ray spectra of AGNs (lAntonuccil 1199311 . 
The torus absorbs the intrinsic X-ray spectra of AGNs (usu¬ 
ally modelled as an absorbed power law) and also scatters 
the X-ray photons to produce a Compton hump at ~20 
keV. The abundant and high quality data from recent X- 
ray satellites enable us to investigate the X-ray properties 
of AGNs in unprecedented detail. The structure of the tori 
in AGNs can be constrained if an X-ray spectral model of 
the tori is specified. Compared with the models for the disk 
geometry, e.g. pexrav an d pexmon (|Magdziarz k Zdziarskil 
119951 : INandra et al.ll2007l l. several more physical and real¬ 
istic models have been recently constructed to model the 
X-ray spectrum from the torus, which adopt a toroidal struc¬ 
ture and self-consistently include fluorescent lines (Ikeda 
et al 2009, hereafter IK09; Murphy & Yaqoob 2009, here¬ 
after MY09; Brightman & Nandra 2011, hereafter BN11). 
These models have been applied to individual AGNs to 
derive the covering factor and column density of the 


sources dSingh et al.ll201 1 : lYaaoobll2012l : lArevalo et al .1120141 : 

iGandhi et al.ll20l4) . With a survey sample, these models are 

helpful for finding Compton thick AGNs and determining 
how the tori evolve with the properties of AGNs, such as 
the correlation between the covering factor and luminos¬ 
ity ( [Brightman fc Nandral 2012^ Brightma n fe Uedal 20121 : 


Bri ghtman et al.1 2014 Buchner et al. 20 14 Ricci et aD 
20l4l . In spite of the success of the application of such mod¬ 


els, the accuracy and validity of such models have not been 
independently verified. The parameters from these models 
will be biased if there are some limitations or errors in these 
models. We recently constructed an X-ray spectral model 
for the tori in AGNs using Geant4, which can handle the 
smooth and clumpy tori by the same code. Thus, if we 
adopt the same assumptions used by MY09, BN 11 and IK09 
(e.g. geometries, cross sections and element abundances), 
our Geant4 code should reproduce their result, in principle. 
In Section 2, we present the quantitative comparison with 
the public models in MY09 and BN11. Because the results 
of IK09 are not public, we just use the general trend pre¬ 
sented in IK09 as a reference. In Section 3, we discuss the 
difference found in the comparison, and we show our main 
conclusions in Section 4. 
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2 SIMULATIONS FOR COMPARISON 

Th e deta ils _of our simulation method are presented in 
iLiu fe Lil J2014h . We have included several physical pro¬ 
cesses in our code, e.g. photoelectric effect, Compton scatter¬ 
ing, Rayleigh scattering, 7 conversion, fluorescent lines and 
Auger process. Similar physical processes have been con¬ 
sidered in the simulations of MY09 and BN11. Thus, it is 
convenient to modify our code to simulate their cases. Our 
code and the simulations by MY09 and BN 11 have included 
multi-scatterings, e.g. Figure 7 in lLiu fe Lil J2014h . We then 
use the same assumptions adopted by BN11 and MY09 re¬ 
spectively to try to reproduce their resulttfj]. The incident 
flux is 5 x 10 8 photons (1-500 keV) for all simulations pre¬ 
sented in this paper. 

2.1 Results under the assumptions of MY09 

The geometry of MY09 is shown in Figure [T] The half¬ 
opening angle of the torus is 60°. 

Using the same geometry, cross sections, element abun¬ 
dances and incident spectra (a single power law with the 
photon index T = 1.8), we can well reproduce MY09’s con- 
tinua, for both Ah = 10 24 and 10 25 cm -2 (Figure [2] where 
di n is the inclination angle relative to an observer). The di¬ 
rect component is not shown for clarity. Since the shape of 
the scattered component is determined by both the scatter¬ 
ing and absorption processes (e.g. Fe K absorption edge at 
7 keV), these results verify the accuracy of MY09’s simula¬ 
tions. 

We have further compared the equivalent width (EW) 
of the F e Kq line of our si mulat ions with that in Fig¬ 
ure 8 of iMurphv fe Yaaoobl (1201111 , which is the erratum 
of MY09. The photon index T is 1.9 in these simulations. 
For Ah = 10 24 and 10 25 cm -2 , the difference of EWs be¬ 
tween our simulations and MY09’s model is smaller than 
1% for most of the directions. The maximum deviation is 
~2% for Ah = 10 24 cm -2 and cos 6% = 0.3 —0.4. This small 
deviation could be due to the different simulation method 
adopted by MY09. The accuracy of this model is sufficient, 
since the statistical error of EW(Fe Ka) in observe d X-rav 
spectra is larger than 10 % in most cases dLiu fe Wandbofol : 
IShu. Yaaoob. fe Wandl20lll ~). 

2.2 Results under the assumptions of BN11 

The geometry of BN 11 is shown in Figure Q5] The half¬ 
opening angle of the torus is fixed at 60°. 

The inner radius of the torus is assumed to be zero in 
BNll’s model. The comparison between the BNll’s model 
and our results is shown in Figure [5] (r=1.8). The total 
spectra are shown, i.e. direct-|- scattered components, since 
BNll’s model only provides the total spectra. The fluores¬ 
cent lines (including the Compton shoulders) in the spectra 
of our simulations are not shown for clarity, because it is 
not easy to directly compare the strength of such a quasi-5 
function in the figure. 

1 The model file of MY09 is downloaded from 
http://mytorus.com/model-files-mytorus-downloads.html 
The model file of BN11 is downloaded from 
http://www.mpe.mpg.de/~mbright/data/torusl006.fits 



Figure 1 . Geometry of MY09. The red and blue arrows indicate 
the possible trajectories of low-energy scattered photons that can 
escape to face-on and edge-on directions, respectively. The yellow 
arrows indicate the trajectories of low-energy scattered photons 
that will be absorbed by the near side of the torus for edge-on 
directions. 



Figure 3. Geometry of BN11. The red and blue arrows indicate 
the possible trajectories of low-energy scattered photons that can 
escape to face-on and edge-on directions, respectively. The yellow 
arrows indicate the trajectories of low-energy scattered photons 
that will be absorbed by the near side of the torus for edge-on 
directions. 

Although the continua of the face-on directions are con¬ 
sistent with each other, the spectra of the edge-on directions, 
especially for low-energy band, are significantly different. It 
seems there is an “additional” component at the low-energy 
band in the simulations of BN11. 

We have further compared EW(Fe Ka) of our simula¬ 
tions with that in Figure 3 of BN11. The photon index T is 
2.0 in these simulations. For the face-on direction (d t0 r = 60° 
and din = 0 — 37°), EW(Fe Ka) of BN 11 is lower than our 
simulation results by 30% and 35% for Ah = 10 24 and 10 25 
cnU 2 , respectively. For the edge-on direction (# t or = 60° and 
d^ = 78 — 90°), their EW(Fe Ka) is higher than our results 
by 60% for Ah = 10 24 cm~ 2 . Due to the large deviation of 
the continuum of the edge-on direction for Ah = 10 25 cm -2 , 
we have not compared the EW(Fe Ka) in this case. 

We explore these discrepancies in detail in the next sec¬ 
tion. 


3 DISCUSSION 

Under the Ah discussed in Figure [2] and [3 the direct com¬ 
ponent is highly absorbed in the low-energy band. Thus, the 
low-energy spectra are dominated by the reflection from the 
visible inner wall of the torus for a given inclination angle, 
which is referred to as the ‘reflection component 2’ in IK09 
(the Figure 2 in IK09 illustrates its geometry and origin). 


















X-ray models for tori in AGNs 3 



Figure 2. Spectra under MY09’s assumptions with JVh — 10 24 (left) and 10 25 cm 2 (right). Only scattered continua are shown for 
clarity. The scattered fluorescent photons are treated as the portion of emission lines and not shown here. 



Figure 4. Spectra under BNll’s assumptions with TVh = 10 24 (left) and 10 25 cm -2 (right). The total spectra (scattered+direct) are 
shown, but the fluorescent lines (including the Compton shoulders) in our simulation are not shown for clarity. The half-opening angle 
is assumed to be 60°. 


The strength of this component depends on the geom¬ 
etry of the torus, i.e. the location and shape of the surface, 
and should also depend on the inclination angles. Therefore, 
to judge whether the results of BNll’s model in Figure[4]are 
reasonable, we investigate the variation of this component 
with inclination angles in different models. 

Because ?'m/r 0 ut = 0.01 adopted in the simulations of 
IK09 is small, IK09’s geometry is actually very similar to 
that of BN11, i.e. the IVh distribution is almost constant for 
different inclination angles. We plot the distribution of IVh 
using equation (3) in IK09 to show this (Figure 0. 

As shown in the Figure 9 in IK09, the ‘reflection com¬ 
ponent 2’ significantly decreases with increasing inclination 
angles, i.e. it is very weak in the edge-on directions. 

A similar trend is also observed in MY09’s model (Fig¬ 


ure 0. though the variation is not as dramatic as that in the 
Figure 9 in IK09. We will explain this later. The spectrum 
at 1 keV with inclination angle 6i n = 65° is higher than that 
with 0i n = 85° by more than one order. 

However, the strength of this low-energy reflection com¬ 
ponent in BNll’s model only weakly depends on the inclina¬ 
tion angles (Figure0. The spectrum at 1 keV with #i n = 65° 
is only higher than that with Sj n = 85° by a factor of two. 

The strength of this low-energy component depends on 
the visibility of the inner wall of the torus at different in¬ 
clination angles. Under the geometry of MY09, since the 
central part of the torus is empty, the whole inner surface of 
the torus is directly illuminated by the central source. As a 
result, a considerable part of the inner wall is visible when 
the inclination angle is slightly larger than the half-opening 
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angle of the torus (the blue arrows in Figure[T|; even for the 
edge-on case, the rim of the inner wall is still visible but the 
majority of reflection in low-energy band is absorbed by the 
near side of the torus (the blue and yellow arrows in Figure 
[lj. This can explain the trend observed in Figure [6] 

To further support the above explanation, in Figure [H] 
we also plot the positions of the scatterings of the observed 
photons in 1-2 keV, i.e. the photons have experienced scat¬ 
terings and finally escaped to the observer. The distribu¬ 
tions of the positions of the scatterings for two directions 
are shown. The scatterings of low-energy photons can only 
occur at the skin of the torus; otherwise, they will be ab¬ 
sorbed in the body of the torus. As the inclination angle 
moves to the edge-on direction, the visible part also moves 
to the rim. If the scattered photons intend to escape from 
the edge-on direction, the scatterings should occur at the 
rim of the torus; otherwise, they will be absorbed by the 
near side of the torus. 

For the BN 11 geometry, the torus extends to the centre 
and the column densities are the same for different incli¬ 
nation angles. As a result, if the low-energy photons are 
scattered and can escape to the observer, they can only be 
scattered very near to the centre; otherwise, they will be 
absorbed before reaching a large radius. Therefore, the re¬ 
flection component is only visible when the inclination is 
slightly larger than the half-opening angle (a few degrees). 
The possible trajectory of the scattered photons is indicated 
by the arrows in Figure[3] For the edge-on case, the scattered 
region is obscured by the near side of the torus. If the scat¬ 
tered photons intend to escape to the edge-on direction, they 
should be scattered at a large radius and follow the way indi¬ 
cated by the blue arrows in Figure^ However, such photons 
should be rare, since most of them will be absorbed before 
reaching the large radius. This is also the reason for the sig¬ 
nificant decrease found in the Figure 9 in IK09. We plot the 
positions of the scatterings in Figure [9j which are indeed 
concentrated in the central part of the torus. The region 
is very small compared with the outer radius (2 pc) of the 
torus. Therefore, this reflection component in the edge-on 
direction should be very weak under the geometry of BN11 
and IK09, as shown by our simulations in Figure 2] 

As mentioned in Section 2.2, the EW(Fe Ka) of BN11 is 
also different from our simulations. For optically thin case, 
EW(Fe Ka) can be analytically calculated, e.g. by equation 
(5) in MY09, and nearly isotropic. Thus, we take this ana¬ 
lytical result as the benchmark test of BNll’s model and our 
simulations. Under the assumption in Figure 3 of BN11, the 
analytical EW(Fe Ka) is 3.7 eV for T = 2.0, #tor = 60°, and 
Ah = 10 22 cm -2 . The EW(Fe Ka) of our simulation under 
the same assumption is well consistent with 3.7 eV (at better 
than 1% level). However, the EW of BNll’s model is about 
4.4 eV for the same case. Since their continua at 6.4 keV are 
consistent with our results at 1% level even for the edge-on 
case of Ah = 10 24 cm -2 , there should be some problems in 
the transportation of Fe Ka photons, which might be related 
to the problem inducing the overestimate of the reflection 
component. 



Figure 5. Ratio between the column density along the line of 
sight and the maximum Ah in IK09 (half-opening angle=60° and 
nn/fout = 0.01). It is almost constant except for the angles near 
the edge of the torus. 



Figure 6. The scattered component is significantly suppressed 
in the edge-on direction under MY09’s geometry. 



Figure 7. The scattered component only weakly depends on 
the inclination angels in BNll’s model. The strength of the low- 
energy component in edge-on direction is much higher than that 
produced by MY09’s model (scattered and direction components 
are added). 
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Figure 8. Positions of the scatterings of the photons (1-2 keV) 
escaped to cos0; n = 0 — 0.2 (top) and cos 0- ul = 0.4 —0.5 (bottom) 
under MY09’s geometry. The 3-D positions are projected onto a 
plane to show the distribution more clearly. The blue line shows 
the boundary of the torus and the central X-ray source is located 
at the origin. 

4 CONCLUSIONS 

With the code constructed using Geant4, we can reproduce 
well the continua and the strength of Fe Ka line of MY09’s 
model. However, the reflection component in the low-energy 
band is much lower than that of BNll’s model. We have 
discussed the origin of this reflection component and shown 
that the scattered region is concentrated in the centre and 
invisible in the edge-on directions under the BNll’s geom¬ 
etry. Therefore, it seems the strength of the reflection com¬ 
ponent is overestimated in BN11 for the edge-on directions. 
The strength of Fe Ka line of BNll’s model is also different 
from our results and the analytical result in the optically 
thin case, which is likely to be due to the problem in the 
transportation of Fe Ka photons. The accuracy of the model 
is crucial to any conclusions from the spectral fitting. 
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